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ABSTRACT 

A formalism for studying the confinement of heavy quarks by considering the 
renormalised quark Dyson-Schwinger equation in the limit txiq — > oo is described. 
We are particularly interested in studying the analytic structure of heavy quark 
propagators in the region of the complex momentum plane close to the bare 
fermion mass pole. 



1. Introduction 

Confinement and chiral symmetry breaking in quantum chromodynamics (QCD) 
must ultimately have a field theoretic cause. They are innately nonperturbative phe- 
nomena, demanding the use of appropriate techniques if they are to be fully understood. 
Lattice gauge theory and the Dyson-Schwinger (DS) equation technique^ represent the 
two most promising lines of attack on QCD from a nonperturbative, field theoretic 
point of view. 

Of these two, lattice gauge theory has the advantage that it is as close as one 
can find to an ab initio treatment of QCD. Approximations are clearly defined and 
in principle controllable. It has the disadvantage that a qualitative understanding of 
physical phenomena tends to be obscured by numerical brute force. Some progress is 
now being made, however, in interpreting lattice gauge fields in terms of monopoles 
and flux tubesB 

Dyson-Schwinger equations provide a complementary approach to QCD. The 
infinite tower of coupled integral equations relating Greens functions is an exact, non- 
perturbative description of QCD. Before any progress can be made, the tower must be 
truncated, generally at the level of quark and gluon propagator equations and bound 
state Bethe-Salpeter and Faddeev equations. Chiral symmetry breaking is signalled 
by the generation of a non-zero scalar part to the quark propagator, while a suffi- 
cient condition for quark confinement is the absence of real, timelike singularities in 
the propagator. These conditions ensue provided the dressed gluon propagator is suf- 
ficiently infrared enhancedB 

Numerical solution of the gluon propagator DS equation is a formidable task, and 
opinions differ about the resultant infrared behaviour of the gluon propagator.™ It has 



been argued however that, if the three gluon and ghost-gluon vertices are modelled in 
a way free of singularities, a strongly infrared enhanced gluon propagator: may result, 
leading to the required chiral symmetry breaking and quark confinement Within the 
DS picture, these two phenomena go hand in hand, both driven by the functional form 
of the quark propagator. 

There is of course no necessary inconsistency between the DS picture and dual 
superconductor flux tube models. A strongly infrared enhanced gluon propagator can 
readily generate a Wilson area law.0 A fortiori, one should not be perturbed that 
the flux tube picture is not immediately manifest in the interpretation of a timelike 
singularity free propagator as a signal of confinement. After all, that Wilson loops 
should satisfy an area law (which we interpret as the existence of a flux tube) is not 
an immediate, or even necessary consequence of the fact that Polyakov loops are zero 
valued (which we interpret as a signal of colour confinement). Obviously, both pictures 
of confinement can happily coexist. 

In what follows we give a very brief descnniion of a formalism designed to extend 
the DS picture to the realm of heavy quarksBlLlJ We shall see that the DS picture of 
quark confinement can be carried over to the extreme heavy quark limit in a natural 
and instructive way. 

2. The Heavy Quark Limit 

Our starting point is the quark DS equation 

Sip)' 1 = t 7 .pA(p 2 ) + B(p 2 ) (1) 
= Z 2 [ 1 -p + m o (A)} + Z 1 J A -0^g 2 D flulfl S(q)T u (q,p), (2) 

with S, D^y and Y ^ the renormalised dressed quark propagator, gluon propagator and 
proper quark-gluon vertex respectively. The renormalisation scale fi is set such that 
S(p)\ p2=/1 2 = l/[«7 ■ V + m R{f J ' 2 )]- We use a Euclidean metric in which timelike vectors 

Satisfy p 2 = -^Minkowski < °> and for Wm ch {l^lu} = 2 <V- 

The heavy quark limit is taken by assuming the expansions 

A(p 2 ,f, 2 ) = 1 + SA( f 2 f + . . . , 5(pV) = m R (ji 2 ) + E B {K, «) + ..., (3) 

where we have defined the momentum variable K and renormalisation point k by 

K = p2 + m «, K =» 2 + m K. (4) 
2imn 2irriR 

The change of independent variable p 2 — > K induced by the transformation Eq. (f|) is 
discussed in detail in reference □ . The main point to note here is that the real timelike 
p 2 axis in the vicinity of the 'would-be bare fermion mass pole', p 2 = —m 2 , becomes 
the imaginary K axis in the new coordinates. The aim of the exercise is to solve for 
the quark propagator in the vicinity of the bare mass pole K — 0. 
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Fig. 1. The modulus |ctq(.K")| of the heavy quark propagator in the complex K plane ob- 
tained by solving the heavy quark DS equation in bare vertex approximation with the gluon 
propagator A(/c 2 ). 

For the purposes of expanding the DS equation in the inverse heavy quark mass 
we introduce the momentum substitutions generic to heavy quark effective theory, 



where v^v^ = 1. This gives K — k ■ v + k 2 /(2im^), and so in the vicinity of k^ = 0, or 
K = 0, the dressed quark propagator has the form, 



The confinement criterion that S(p) should be free from timelike singularities on the 
negative real p 2 axis translates in the heavy quark case to a requirement that gq should 
be free from singularities on the imaginary K axis. 

In reference the heavy quark propagator is solved for using the bare vertex 
or "rainbow" approximation r^(g,p) = 7^, and Landau gauge smeared Frank and 



P^ = irriRV^ + k i 



(5) 




(6) 



RobertsS model gluon propagator D^ u (k) = (8^ — k fl k u /k 2 )A(k 2 ), where 

a<* 2 ) = m 4 ^-"" + — • ro 

Here d = 12/(33 — 2Nf), Nf = 3 is the number of light quark flavours, rrit = 0.69 GeV 
is a parameter fitted to a range of calculated pion observables, and the Gaussian width 
a is chosen to be 0.5643 (GeV) 2 . 

In Fig. |l] is plotted the modulus \ctq(K)\ of the heavy quark propagator obtained 
by solving the DS equation in the vicinity of K = 0, with a renormalised heavy quark 
mass of rrtR = 5 GeV. We see that the bare propagator pole at K = has moved 
well away from the imaginary K axis, leaving the heavy quark propagator analytic 
over a large region in the vicinity of K = 0. The removal of propagator poles from the 
imaginary K axis signals heavy quark confinement. As noted in reference □ , a large shift 
of the pole portends well for the possibility of finding solutions to the Bethe-Salpeter 
equation for Q-q mesonsB Solution of the heavy meson Bethe-Salpeter equation, of 
which accurate determination of heavy quark propagators is a necessary first step, 
should prove helpful in impulse approximation studies of semileptonic decays.^ 2 
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